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Section  1 


INTRODUCTION 


1 . 1  UNEAR  CHAIR  REQUIREMENTS 

Project  Linear  Chair  is  examining  methods  of  measuring  the  near 
field  and  far  field  magnetic,  electric,  and  electro-magnetic  signatures  of 
various  underwater  platforms.  Initial  hydro-mechanical  mooring  approaches 

which  have  been  proposed  to  position  the  sensor  strings  ari' depicted  In- - 

Plgure-i-l  for  the  near  field  system  and  Ptgare'T'^^^for^e  far  field  system . 

A  stringent  requirement  of  any  proposed  mooring  system  is  the 
movement  of  the  attached  sensors  in  response  to  ocean  currents.  A  sensor 
watch  circle  of  3  ft.  radius  is  a  desired  goal  to  achieve. 

1 . 2  OVERVIEW  OF  REPORT 

The  following  report  contains  analysis  methods  for  determining 
the  static  deflections  that  could  be  expected  from  the  mooring  configurations 
depicted  in  Figures  1-1  and  1-2.  The  analysis  methods  have  been  exercised 
for  various  parametric  combinations  of  system  components  to  obtain  an 
initial  range  of  excursions  that  could  be  expected  under  different  current 
conditions.  It  is  emphasized  that  these  initial  results  are  primarily  presented 
for  illustrative  purposes  since  gross  assumptions  were  made  concerning  the 
size  and  weight  of  system  components.  However,  as  the  specific  designs  are 
evolved,  the  analysis  techniques  presented  herein  can  serve  as  a  valuable 
design  tool  to  assess  the  hydro-mechanical  performance  to  be  expected. 


Section  2  addresses  the  near  field  array  analysis  and  Section  3 
addresses  the  far  field  (goal  post)  array  analysis.  Appendices  have  been 
used  liberally  to  present  detailed  mathematical  formulations  and  descriptions 
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Candidate  Mooring  Configurations  for  the  Near-Field  Array  System 


Section  2 


NEAR- HELD  ARRAY  ANALYSIS 
2 . 1  CONnCURATION  DISCUSSION 


The  proposed  near  field  array  consists  of  two  independent  sensor 
strings  separated  by  approximately  180  ft.  Table  2-1  contains  the  pertinent 
parameters  examined  in  this  analysis. 


Table  2-1 .  Configuration  Parameters  for  Near- Field 
Array  Analysis 


Water  Depth 
Cable  Parameters 
Construction 
Diameters 
Overall  length 
Buoyed  length 
Wt/ft,  (water) 
Buoy  Parameters 
Material 
Net  Buoyancies 
Diameters 

Shape 

Drag  Coefficient 
Currents 
Type 

Magnitude 


3000  ft. 

Kevlar 

0.3,  0.5,  and  1.0  inches 
2600  ft. 

1000  ft. 

0  (neutrally  buoyant) 

Aluminum  6061-T6 

500,  1000,  3000,  5000  lbs.  with 

of  3.02,  3.80,  5.49,  6.51  ft. 
respectfully 

Spherical 

0.5  based  on  frontal  area 

Uniform  for  entire  depth 

0.05,  0.1,  0.3,  0.5,  and  1.0  kts 


The  buoy  parameters  were  calculated  based  on  a  3000  ft.  design 
depth  and  a  safety  factor  of  3.  Detailed  calculations  of  the  buoy  parameters 
are  contained  in  Appendix  A. 


The  configurations  analyzed  were  of  two  types  and  are  depicted  in 

Figures  2-1  and  2-2.  In  the  two-buoy  configurations  the  size  of  the  buoy  at 

the  top  and  intermediate  positions  were  sequentially  changed.  The  resultant 
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Top  buoys  with  net  buoyancies 
of  500,  1000,  3000  and  5000 
pounds 

Diameters  of  3.02,  3.80,  5.49 
and  6.51  feet 


Intermediate  buoy  with 
Net  Buoyancies  of  500, 
1000,  3000  and  5000  pounds 
diameters  of  3.02,  3.80, 
5.49  and  6.51  feet 


mT7V 


Uniform 
Currents 
of  0.05, 
0.1,  0.3, 

0.5  and 

1 .0  knots 

Water 
Depth  of 

3000  feet 


Figure  2-1 .  Two-Buoy  Near-Field  Array  Analysis  Diagram 


All  buoys  have  500  pounds 
Net  Buoyancy  with  3.02  foot 
diameter 


0,5  inch  Smooth  Neutrally 
Buoyant  Kevlar  Cable 


1000  ft 


1600  ft 


mrrrm 


Uniform 
Currents 
of  0.5 
and  0.3 
knots 


Water 

Depth 

of 

3000  feet 


Figure  2-2.  Distributed  Buoy  Near-Field  Array  Analysis  Diagram 


configurations  were  then  analyzed  in  uniform  currents  of  0.05,  0.1,  0.3, 


0.5,  and  1.0  knots.  In  the  distributed  buoy  configuration  11  500-lb.  buoys 
were  uniformly  spaced  over  the  1000-ft.  buoyed  length.  This  configuration 
was  analyzed  in  currents  of  0.3  and  0.5  knots. 

2 . 2  RESULTS  OF  TWO-BUOY  CONHOURATION 

Tables  2-2  thru  2-5  contain  the  computed  watch  circle  radii  for  the 
top  and  intermediate  buoy  under  different  configurations  and  currents.  The 
tables  are  for  net  buoyancies  of  500,  1000,  3000,  and  5000  lb  at  the  top, 
respectfully. 

The  combinations  that  satisfy  the  initial  watch  circle  radius  of 
3  ft.  for  the  top  buoy  are  delineated  by  the  red  line  in  the  tables. 

2.3  RESULTS  OF  DISTRIBUTED  BUOY  CONFIGURATION 

In  addition  to  examining  a  two  buoy  approach,  uniformly  distributed 
loading  of  the  top  1000  feet  of  the  cable  was  also  investigated.  The  following 
parameters  were  used: 

Cable 

Diameter  0.5  inch 

Wt/ft.  (water)  0  (neutrally  buoyant) 

Buoy 

Net  Buoyancy  500  lbs. 

Diameter  3.02  ft. 

Number  11  buoys  uniformly  spaced  over 

top  1000  ft.  scope. 

Current 

Type  Uniform 

Magnitude  0.3,  0.5  knots. 

The  Intent  of  examining  the  uniformly  spaced  buoyancy  was  to  obtain 

tradeoffs  In  watch  circle  performance  for  distributed  versus  concentrated 

buoyancy  elements  in  a  uniform  current.  The  cases  compared  are  depicted 

graphically  in  Figure  2-3.  The  composite  net  buoyancy  in  each  system 
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821.59  6J7.96  681.92  498.29  438.50  254.87  358.44  171.81 


e  2-4.  Watch  Circle  Radii  with  3000  Pounds  Net  Buoyancy  at  the  Top  for  Various 
Intermediate  Position  Buoyancies  and  Uniform  Currents 


compared  totaled  5500  lb.  As  observed  in  the  data,  the  two-buoy  system  with 
the  leirger  (5000-lb)  buoy  at  the  top  results  in  the  better  watch  circle  performance 
when  the  configuration  is  exposed  to  a  uniform  cunrent.  This  result  is  generally 
true  for  uniform  currents;  however,  for  sharp  current  profiles,  a  distributed  buoy¬ 
ancy  approach  is  generally  better  because  the  buoys  are  placed  in  lower  currents 
and  thus  result  in  reduced  system  drag.  When  current  data  from  specific  sites 
are  known,  more  meaningful  tradeoffs  can  be  made. 
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FAR-FIELD  ARRAY  ANALYSIS  (GOAL  POST) 

3 . 1  OVERVIEW 

The  goal  post  array,  in  the  analytical  sense,  represents  a  system 
which  must  be  analyzed  in  an  iterative  fashion.  Figure  3-1  depicts  the 
configuration  which  has  been  examined.  The  intent  has  been  to  obtain  a 
technique  that  can  determine  the  absolute  current  response  of  various 
system  components  relative  to  a  zero  current  condition .  For  the  material 
presented  herein  the  current  has  been  assumed  to  be  uniform  and  planar; 
however,  the  technique  can  easily  be  extended  to  cover  currents  at  various 
azimuths . 

3 . 2  ANALYSIS  APPROACH 

The  basic  steps  required  to  analyze  the  configuration  in  Figure  3-1 


is  as  follows: 


1 .  Choose  basic  system  components  and  determine  the 
zero  current  equilibrium  configuration. 

2 .  Analyze  the  response  of  the  two  vertical  strings  relative 
to  their  connection  point  on  the  array.  Since  the  vertical 
strings  have  a  free  end  their  relative  motion  will  be 
independent  of  the  movement  of  the  trapezoid.  This  is 
strictly  true  for  uniform  currents  and  may  or  may  not  be 
true  for  current  profiles.  The  latter  case  is  dependent  on 
the  vertical  response  of  the  trapezoid  in  that  it  may  position 
the  vertical  string  in  a  different  current  regime. 

3 .  The  analysis  of  (1)  will  output  a  tension  vector  at  the 
trapezoid  apex  which  will  be  used  in  the  trapezoidal 
configuration  analysis . 

4.  Analyze  the  trapezoidal  array  using  Newtons  method  (described 
in  Appendix  B)  to  determine  the  static  response  configuration. 


ORIGIN  OF  COORDINATE  SYSTEM 


5.  Combine  the  geometry  response  of  the  vertical  strings 
and  the  trapezoid  to  obtain  the  overall  response  of  the 
vertical  string  top  buoy  which  is  to  remain  as  stationary 
as  possible. 

The  following  sections  discuss  the  specifies  of  each  step  discussed 

above. 

3 . 3  BASIC  CONFIGURATION 

Figure  3-2  depicts  the  equilibrium  diagram  for  the  zero  current  case 
that  has  been  assumed  for  illustrative  purposes.  The  forces  and  angles  result 
from  straightforward  calculations.  The  trapezoidal  cables  are  straight  lines 
since  it  is  assumed  that  they  are  neutrally  buoyant  with  no  hydrodynamic 
forces  present  from  currents. 

3.4  ANALYSIS  OF  VERTICAL  STRINGS 

The  vertical  strings  were  analyzed  parametrically  for  various  buoys , 
cable  diameters  and  current  velocities.  The  resultant  watch  circle  data  represents 
a  "best  case"  in  total  watch  circle  response,  because  it  assumes  the  trapezoid 
to  be  stationary.  Figure  3-3  depicts  the  parameters  assumed  for  the  analysis. 

The  technique  is  similar  to  that  used  on  the  near  field  array  and  discussed  in 
Section  2. 

Table  3-1  contains  the  analytical  results  in  terms  of  watch  circle 
radius  and  vertical  displacement  of  the  top  buoy  relative  to  the  connection  point 
at  the  apex  of  the  trapezoid.  Figure  3-4  graphically  depicts  the  same  results. 

Note  that  the  design  goal  of  3-foot  watch  circle  radius  is  highlighted  and  eliminates 
the  majority  of  configurations  examined. 

For  purposes  of  the  goal  post  array  system  analysis  the  solution  using  a 
5000  lb.  buoy  and  a  1 .0  inch  diameter  Kevlar  cable  in  a  1-knot  current  was  used. 
From  the  analysis  of  this  case  the  vertical  string  tension  at  the  trapezoid  was 
5997.42  lbs.  tension  at  an  angle  of  95.33°  degrees. 


Equilibrium  Diagram  for  the  Trapezoidal  Array 
With  Zero  Current 


PRELIMINARY  WATCH  CIRCLE  TRADE-OFFS 

ft 

OF  VERTICAL  SENSOR  STRINGS 


Figure  3-4.  Preliminary  Watch  Circle  Trade-offs  of  Vertical  Sensor  Strings 
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Geometric  Response  for  Vertical  Strings  for  Various  Combinations 
of  Buoyancy,  Cable  Diameter,  and  Uniform  Current  (cont.) 
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3 . 5  DETERMINE  RESULTANT  FORCE  ON  TRAPEZOID 

In  addition  to  the  vertical  string  tension  vector  it  was  assumed 
that  a  1000  lb.  buoy  was  located  at  each  trapezoid  apex.  Figure  3-5  depicts 
the  vector  force  diagram  that  computes  the  resultant  force  acting  on  the 
trapezoid. 

3 . 6  ANALYSIS  OF  THE  TRAPEZOID 

The  application  of  Newton's  method  in  calculating  the  catenary  of 
the  goal  post  configuration  current  response  consists  of  utilizing  the  known 
zero  current  geometry  of  the  goal  post  array  and  utilizing  the  catenary 
computer  program  to  successively  calculate  trial  solutions.  The  process  is 
started  by  (1)  assuming  a  new  tension  vector  at  the  base  of  the  trapezoid; 

(2)  assuming  a  perturbed  set  of  tension  vectors  about  that  chosen  in  (1); 

(3)  calculating  the  resultant  catenary  of  the  trapezoid;  (4)  comparing  the 
catenary  results  to  the  compatibility  equations;  (5)  predicting  the  next 
initial  condition  from  the  error  derivatives;  and  (6)  computing  the  system 
watch  circles.  The  following  example  Illustrates  this  procedure  whose 
theoretical  approach  is  detailed  in  appendix  6. 

3.6.1  Assume  a  New  Tension  Vector 

Figure  3-6  depicts  an  intuitive  response  geometry  to  the  trapezoid 
from  an  imposed  current.  The  tension  will  increase  and  the  angle  will  increase 
at  Point  A.  Therefore,  an  assumed  solution  is; 

Tj  =  16,200  lbs  (Tq  was  16,163.57) 

0^  =  159.0  degrees  was  158.19) 

3.6.2  Assume  a  Perturbed  Set  of  Initial  Ctonditions 

The  purpose  cf  this  is  to  find  the  effect  of  small  perturbations  of  the 


assumed  solution  on  the  compatibility  relationships  in  order  to  form  error 
derivatives. 


Figure  3-6.  Intuitive  Response  Geometry  to  a  Flaneur  Current  for  the  Goal  Post  Array 


For  purposes  of  this  example  the  two  sets  are: 


+AT  =  16,200  +  5  =  16,205  lbs 


=  159  degrees 


Tj  =  Tj  =  16,200  lbs 

=  159  +  (0.1)  =  159.1  degrees 


3.6.3  Calculate  the  Resultant  Catenary  of  the  Trapezoid 
Using  the  three  sets  of  initial  conditions: 

Ti , 

Tj  +AT,  jar^ 

T,  0^  +^0 

The  resultant  geometry  of  the  end  of  the  trapezoid  relative  to  the 
known  coordinates  of  the  end  point  (0,0)  are:  (See  Figure  3-7) 


Initial  Condition 
Tj  =  16,200  lbs 

fifj  =  159° 


23.24 


-137.78 


Tj  =  16,205  lbs 
0^  =  159° 


22.90 


-136.95 


Tj  =  16,200  lbs 
0^  *  159.1° 


22.94 


-149.29 


3.6.4  Compute  Errors  from  the  Compatibility  Relationships 

From  the  known  relative  geometry  of  the  anchor  points  of  the  trapezoid 
the  following  conditions  must  hold  (refer  to  Figure  3-7): 

f  =  X  =  0 


g  =  y  =  0 


3-13 


v  :4W’ 


Trapezoidal  Array  Nomenclature 


For  the  cases  above  the  error  derivatives  are:  (see  Appendix  B  for  derivatives) 


3.6.5 


gf  _  22.90  -  23.24 
dT  ~  16,205  -  16,200 


=  -0.068 


df  ^  22.94  -  23.24  _  , 
a0  159.1  -  159 

fa  ^  -136.95  -  (-137.78)  _  «  ififi 
el  16,205  -  16,200  "  * 

§3.  _  -149.29  -  (-137.78)  _ 

30'  159.1  -  159  -’115.10 

Predict  the  Next  Initial  Condition  for  T  and  0  in  Order  for  f  = 
The  next  set  of  conditions  are: 


=  Tj  +dT 

j  +  d0 


where  =  16,200  lbs 


dT  =  - 


—30 

3f  3q 
3T  d0 


^  30 
3Q  3l 
5T  d0 


dgr  =  - 


5  .  (  22 
at  ^  ST 

3i  dQ  _  3Q  3f 

dT  00  ~  3T  30 


Using  the  numbers  calculated  in  section  3.6.4  the  new  assumptions 
for  T  and  0  are: 


T^  =  16,200 


‘(23.; 


24)  (-115.10)  -  (-137.78)  ( 
.068)  (-115.10)  -  (0.166)  ( 


(-3.00)' 
-3.00)  _ 


=  16,200  +  370.97  =  16,570.97 
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02  =  159  + 


of  =  I«;q  +  (“137.78)  (-0.068)  -  (23.24)  (0.166) 

^2  (-0.068)  (-115.1)  -  (0.166)  (-3.00) 

=  159  -  0.66  =  158.34 

A  new  set  of  perturbed  conditions  about  ,  0^  were  computed  and 
the  process  continued.  After  a  total  of  three  iterations ,  the  following  results 
were  obtained  for  a  1  knot  uniform  current: 


xm 


Y  (ft) 


Tension  (lb) 


below: 


0 

7100.00 

0 

16,597.46 

158.38 

3500 

3847.55 

1292.85 

16,561.21 

157.74 

3500 

3847.55 

1292.85 

14,772.90 

178.85 

4100 

3247.67 

1304.90 

14,770.65 

178.85 

4100 

3247.67 

1304.90 

15,303.77 

158.21 

7600 

0.06 

-0.02 

15,267.23 

157.51 

above  data  can  be  compared  to  the  zero  current  case  presente 

Scope  (ft) 

X  (ft) 

Y  (ft) 

Tension  (lbs) 

0  (degrees 

0 

7100.00 

0 

16163.6 

158.19 

3500 

3850.00 

1300.00 

16163.6 

158.19 

3500 

3850.00 

1300.00 

15006.6 

180.00 

4100 

3250.00 

1300.00 

15006.6 

180.00 

4100 

3250.00 

1300.00 

16163.6 

158.19 

7600 

0 

0 

16163.6 

158.19 

3.6.6  Compute  the  System  Watch  Circles 

The  resultant  watch  circle  of  the  various  parts  of  the  goal  post 
array  system  are  computed  below: 

•  Displacement  of  Buoy  "A"  in  a  1  knot  current  relative  to  a 
zero  current  condition: 

X  displacement:  3847.55  -  3850  =  -2.45  ft 
y  displacement:  1292.85  -  1300  =  -7.15  ft 
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•  Displacement  of  Buoy  "B"  in  a  1  knot  current  relative  to  a 
zero  current  condition: 

X  displacement:  3247.67  -  3250  =  -2.33  ft 
y  displacement:  1304.90  -  1300  =  4.90  ft 

•  Displacement  of  Vertical  string  buoy  "C"  relative  to  buoy 
"A"  and  buoy  ”D"  relative  to  buoy  "B": 

X  displacement  =  75.69  ft 
y  displacement  =  -2.62  ft 

•  Total  displacement  of  buoy  "C: 

X  displacement  =  -75.69  -  2.45  =-78.14 
y  displacement  *  -  2.62  -  7.15  =  -9.77  (dip) 

•  Total  displacement  of  buoy  ”D": 

X  displacement  =  -75.69  -  2.33  =-78.02 
y  displacement  =  -2.62  +  4.90  =  2.28  (rise) 

3 . 7  OBSERVATIONS  OP  FAR- FIELD  GOAL  POST  ANALYSIS 

The  analysis  presented  in  the  preceding  sections  have  demonstrated 
a  method  of  assessing  the  array  response  to  the  environmental  effects  of 
current.  It  is  apparent  from  the  analysis  of  the  vertical  strings  (Figure  3-4) 
that  watch  circles  in  the  order  of  3  ft.  can  only  be  maintained  under  very  low 
current  conditions  (.4  knots  or  less  depending  on  buoyancy  and  size  of  cable). 
The  analysis  model  presented  herein  has  been  intended  to  demonstrate  an 
approach  for  estimating  array  response  to  ocean  currents.  Definition  of  system 
components  and  in- situ  current  profiles  are  needed  to  provide  more  accurate 
results;  however,  the  model  should  be  used  in  a  parametric  sense  to  assist  in 
arriving  at  a  final  design.  Finally  the  model  can  be  easily  extended  to  examine 
non-planar  currents  for  a  total  excursion  analysis  capability. 


From  Table  3-1 
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CALCULATIONS 
OF  SPHERICAL  BUOY 
PARAMETERS 

The  buoyant  member  was  assumed  to  be  a  spherical  buoy  made  of 
6061-T6  aluminum  with  a  3000  foot  working  depth  and  a  design  safety  factor 
of  3. 

Computations  were  made  to  determine  the  wall  thickness  versus  the 
inner  and  outer  radius  of  the  spherical  buoy.  An  incremental  element  of  stres 
was  defined  as 

<7  =  qR2/2t 

where  q  =  the  external  pressure 

R2  =  the  inner  wall  radius  of  the  buoy 
t  =  the  buoy  wall  thickness 

and  q  =  P  g  h 

=  (1.99  sluq/foot^)  (32.2  ft/sec^  )  (3000  feet) 

144  inch^/foot^ 

=  1335  p.s.i. 

The  maximum  allowable  stress  is  given  by 

a  max  =  Compressive  Yield  Strength  of  Material 

Safety  Factor 


For  6061-T6  aluminum  the  compressive  yield  strength  is  35,000  psi.  Therefore 

35000 

crmax  =  =  11667  psi 


Thus 


Since 


Then 


a  max  =  ^  *• 

R^/t  =  2(11667)/1135  =  17.5 


=  R  -  t/2 


where  R  =  the  outer  radius  of 
the  buoy 


R  -  t/2 


=  17.5 


R  =  17, 5 1  +  t/2 
R/  =  18 

Calculations  were  then  made  of  the  net  buoyancy  versus  radius.  Since  the  net 
buoyancy. 


Then 


Bnet  =  Bgross  -  Weight  of  the  buoy 

Bnet  =  V3  ^  R^  7  water  -  |^3  R^  7al  -  tt  (R-t^)  7al 

=  Vs  .R^(  7  water  -  7al)  -  (R-t)^  7al| 

Since  R/t  =  18 

Bnet  =  V3  tt  R^(7water  -  7al)  -  (R-Vl8)^  7al 


=  /3  ^[R  (  7 water  -  7al)  -  1.176  R  7al 

-  Vs  ^  R^(  7water  -  0.1761  7al) 

3 

Since  7water  =  64.4  Ib/ft 
7al  =169.3  Ib/ft^ 


Bnet  =  144.8  R' 


where  R  is  in  feet 


Subsequently  Table  A-1  presents  the  buoyancies  examined  versus  the 
outside  diameter  of  the  aluminum  spheres  required  to  produce  such  buoyancy. 


Table  A-1  Net  Buoyancy  and  Buoy  Diameter 


B(net) 

(lbs) 

Diameter 

(feet) 

50 

1.40 

100 

1.76 

300 

2.55 

500 

3.02 

1000 

3.80 

3000 

5.49 

5000 

6.51 

The  semifixed  array  was  evaluated  in  uniform  currents  of  0.05,  0.1, 
0.3,  0.5,  and  1.0  Knots.  Consequently  the  drag,  D,  caused  by  the  buoy  was 
computed  by 

D  =  (l/2pV^) 

where  =  drag  coefficient 

=  .5  for  10^<  Re  <  3  X  10^ 

A  =  spherical  area 
-  TrD^4 

P  =  Density  of  Water 
=  1.99  slugs 

=  Velocity  in  Knots^^l.689  ft/sec  j  ^ 
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APPUCATION  OF  NEWTON'S  METHOD 
IN  CALCULATING  THE  EQUILIBRUM 
CONHCURATION  OF  A  TWO-LEGGED 
BUOY  MOORING  WITH  AN  IN  PLANE  CURRENT 

Newton's  method  is  basically  a  numerical  iterative  means  for  finding 
a  zero  of  a  function  ftx).  This  method  can  be  extended  to  mathematical  systems 
of  n  equations  in  n  unknowns . 

For  application  to  the  trapezoided  array  reference  is  made  to  Figure 
B-1 .  In  order  to  calculate  an  assumed  configuration  due  to  a  current  a  new 
tension  vector  at  point  0)  is  assumed  and  the  resultant  location  of  the 
end  point  (f,  g)  determines  how  accurate  your  initial  assumption  was.  In 
order  for  your  assumption  to  be  correct: 

f  =  g  =  0  (B-1) 

The  total  differential  of  f  and  g  can  be  expressed  as: 

“  If  It  ■  4°  -  «i 

In  order  for  and  to  be  equal  to  zero,  the  following  expression 
in  matrix  form  must  hold: 


5T 

a9 

dT 

i3. 

aT 

is. 

ae 

d6 

or. 


B-1 


CURRENT  CONDITION 


Figure  B-1 .  Trapezoidal  Array  Nomenclature 


d0  =  - 


il 

^ _ L 

gf  gq 

gT  ge 


f  ^ 

gT 

gf  gq 

g0  aT 


(B-5) 


dT  =  - 


I  S3.  .  Si 

gQ  ^  ge 

^52.  gf  gg 
gT  ge  “  ge  gT 


(B-6) 


Therefore  the  next  choice  for  T_  and  ©2  order  to  make  f.  =  g,  =  0  are: 


T2  =  Ti  +  dT  (B-7) 

02  =  e^  +  de  (B-6) 

where  dT  and  de  are  defined  in  equations  B-6  and  B-5  respectfully. 

The  error  derivatives  and  are  evaluated  by  obtaining 

values  of  f  and  g  for  very  small  variations  of  the  initial  assumption  Tj  ,  0^  , 


that  is: 


gf  -  ^(T^.  9^) 

gg  “  AT 

ga  .  ^(Ti  e,)  -  9(11,  9^) 
gT  '  AT 

91  *(Ti.  -  ^CTi.  V 


ia  - 
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9(Tj,  0^  +  A0)  -  9(Tj^9  0^) 


(B-7) 


(B-8) 


(B-9) 


(B-10) 


The  solution  will  converge  quadratically  and  as  in  the  example  presented 
in  Section  3,  takes  about  4  iterations  in  order  to  obtain  f  =  g  =  0. 


